PT-symmetric double-well potentials revisited: bifurcations, stability and dynamics 
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In this work we analyze PT-symmetric double-well potentials based on a two-mode picture. We 
reduce the problem into a PT-symmetric dimer and illustrate that the latter has effectively two 
fundamental bifurcations, a pitchfork (symmetry-breaking bifurcation) and a saddle-center one, 
which is the nonlinear analog of the PT-phase-transition. It is shown that the symmetry breaking 
leads to ghost states (amounting to growth or decay); although these states are not true solutions of 
the original continuum problem, the system's dynamics closely follows them, at least in its metastable 
evolution. Past the second bifurcation, there are no longer states of the original continuum system. 
Nevertheless, the solutions can be analytically continued to yield a new pair of branches, which is 
also identified and dynamically examined. Our explicit analytical results for the dimer are directly 
compared to the full continuum problem, yielding a good agreement. 



I. INTRODUCTION 

PT-symmetric quantum systems [H 0] have emerged 
as an intriguing complex generalization of conventional 
quantum mechanics and have been a focus for numer- 
ous investigations at the interface between theoretical 
physics and applied mathematics. The key premise is 
that fundamental physical symmetries such as parity P 
and time reversal T may be sufficient (in suitable para- 
metric regimes) to ensure that the eigenvalues of the 
Hamiltonian are real. Thus, PT-symmetric Hamiltoni- 
ans provide an alternative to the standard postulate that 
the Hamiltonian operator be Dirac Hermitian (invariant 
under matrix transposition and complex conjugation *). 
In the context of Schrodinger Hamiltonians with a com- 
plex potential V(x), the constraint of PT symmetry re- 
quires that the potential satisfy V(x) — V*(— x)\ that 
is, V(x) has a symmetric real part and an antisymmetric 
imaginary part. 

On the other hand, nonlinear Schrodinger (NLS) 
equations incorporating double-well potentials have re- 
ceived attention due to applications in atomic and op- 
tical physics. Such potentials can easily be realized in 
the context of atomic Bosc-Einstcin condensates (BECs) 
through the combination of a parabolic (harmonic) trap 
with a periodic potential. Experiments have observed 
fundamental effects, including tunneling and Joseph- 
son oscillations for small atom numbers, macroscopic 
quantum self-trapped states for large atom numbers 
|3j , and nonlinearity- induced symmetry-breaking dynam- 
ical instabilities Theoretical studies accompanying 
these developments have examined finite-mode reduc- 
tions, analysis of the bifurcations and their dynamical 
implications [H-EIJ, as well as quantum effects [p| and 
nonlinear variants of the potentials [l4j|. A similar phe- 
nomenology has also been found in nonlinear optical set- 
tings, with results for the formation of asymmetric states 
in dual-core fibers [IH, self-guided laser beams in Kerr 
media [lfjj . and optically-induced dual-core waveguide 



structures in photorefractive crystals [171 ] . 

Recently, double-well potentials in the context of PT- 
symmetric nonlinear systems have received consider- 
able attention. This is due to the pioneering work of 
Christodoulides and collaborators, who proposed that 
nonlinear optics presents a fertile ground for the exper- 
imental realization of PT-symmetric systems. The first 
realization of PT-symmctry in a waveguide coupler arose 
in the so-called passive-PT setting in which two waveg- 
uides, one with loss and the other without loss, were used 
[l8| . A similar proposal for the existence of a leaking 
dimer (in the presence of nonlinearity) was formulated in 
the atomic setting of Bose-Hubbard models [l|| . Subse- 
quently, an optical waveguide system with both gain and 
loss was studied and the role of the nonlinearity in the 
dynamics was explored [2fj| . Further experimental inves- 
tigations were concerned with electrical analogs of the 
system [2l[. Theoretical investigations have rapidly fol- 
lowed by examining such dimer-type settings j22H30{ and 
generalizations thereof, including ones where the gain- 
loss contributions appear in a balanced form in front of 
the nonlinear term [3lM33| . 

This paper revisits this theme of PT-symmetric 
double-well potentials. Our motivation is to unify the 
above studies with an important recent contribution, 
namely, Ref. [34| . To be specific, in the early works on 
NLS models with double- well potentials d [ill, El HI , 
and also in the PT-symmetric dimer (28[ , the symmetry- 
breaking bifurcation was identified, but the asymmetric 
states that normally accompany the bifurcation [El [Til Il2j 
could not be identified. The simpler system of a PT- 
symmetric pair of 5-function potentials, where the so- 
lution can be obtained by means of a five-dimensional 
numerical root search, was studied in Ref. [33 ]. There, it 
was found that the bifurcation results in the emergence 
of what we call a ghost state, namely, a solution of the 
steady-state problem with a complex nonlinear eigenvalue 
parameter (complex propagation constant in optics or 
complex chemical potential in BECs). [These states were 
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discussed in a 2006 paper for the case of the (discrete) 
leaking dimer Another observation of Ref. [34[ con- 
cerned the possibility of performing an analytic continua- 
tion of the symmetric and antisymmetric solutions of the 
original double-well problem past the point of their VT 
phase transition; that is, past the critical point of the 
saddle-center bifurcation where they disappear simulta- 
neously. This is in the broader spirit of PT-symmetric 
quantum mechanics @] . 

In this paper we study an NLS model with a double- 
well potential, which has an even real part and a 'PT- 
symmetric imaginary part. We combine the above ideas, 
interweaving analytics based on the two-mode analysis of 
the simple dimer model and numerics in the framework 
of the partial differential equation (PDE). 

This paper is organized as follows. We present the 
two-mode analysis of the discrete nonlinear Schrodingcr 
equation (DNLS) dimer in Sec. II, where we examine the 
symmetric/antisymmetric and asymmetric states and the 
potential for nonlinearity-induced symmetry-breaking bi- 
furcations. We show how the 'PT-symmetric variant of 
the dimer problem adjusts these features. The modified 
critical point of the pitchfork bifurcation is identified, 
and the feature of a saddle-center bifurcation, which is a 
nonlinear analog of the VT phase transition, is discussed 
in the spirit of some of the earlier work on this theme. 
The key by-products of these bifurcations, namely, the 
ghost states emanating from the pitchfork and the an- 
alytically continued states past the nonlinear analog of 
the VT phase transition, are found in analytical form. 
Armed with these analytical results, we study the full 
PDE problem in Sec. Ill, where the aim is to reduce 
the PDE, through suitable and fully quantified approxi- 
mations and transformations (so that one can translate 
the dimer results to the PDE and vice- versa), to a PT- 
symmetric dimer. In Sec. IV we compare the analyti- 
cal results concerning the dimer problem with numeri- 
cal ones obtained in the framework of the original NLS 
model. In addition to quantifying the bifurcations and 
identifying their "daughter" states, we examine the dy- 
namics of the various unstable states within the system. 
Finally, in Sec. V we briefly summarize our findings and 
present some possibilities for further study. 

II. DNLS DIMER 

A. Classical DNLS dimer 

We begin our exposition by revisiting a simpler initial 
problem that has been long studied (see, e.g., Ref. (HI), 
namely the DNLS dimer of the form 

iii\ = — ku2 — |ui| 2 Mi, 

iii 2 = —kui - \u 2 \ 2 u 2 , (1) 



tance in optics). As is customary, we seek stationary so- 
lutions of the form u\ = e~>q>(iEt)a and u 2 = exp(iEt)b, 
where the complex amplitudes a and b and the propaga- 
tion constant E are determined by the equations 

Ea = kb+ \a\ 2 a, Eb = ka+\b\ 2 b. (2) 

Using a polar representation of the two "sites" , namely 
a = Ae l0a , b = Be l6b , we immediately obtain sin (9 = 0, 
where 9 = 9f, — 9 a , and hence the resulting equations for 
the amplitudes are 

EA = ±kB + A 3 , EB = ±kA + B 3 . (3) 

The simpler solutions of the above system are symmetric 
ones of amplitudes A 2 = B 2 = E =p k with the upper 
(lower) sign corresponding to the in-phase (out-of-phase) 
profile. In addition, there exist asymmetric solutions of 
amplitudes A =/= B, which can be determined by the fol- 
lowing equations derived by eliminating E from ([5]): 

AB = ±k, A 2 + B 2 = E. (4) 

The amplitudes A 2 = (E ± VE 2 ~Ik 2 )/2, B = k/A ex- 
ist only for E 2 > \k 2 and thus for E > 2k or E < -2k. 
These asymmetric solutions coincide with the symmet- 
ric ones at E 2 — 4k 2 , and hence a pitchfork symmetry- 
breaking bifurcation is responsible for their emergence 
in the nonlinear problem. As is known from previous 
works [H, [ll|, Ell ■> this bifurcation arises for the focusing 
(attractive) nonlinearity case k > at the point E = 2k, 
and does so from the symmetric branch. However, in the 
defocusing (repulsive) case k < the bifurcation arises 
when E = —2k and the asymmetric solution emerges 
from the antisymmetric branch. This bifurcation picture 
is also complemented by the stability eigenvalues for the 
focusing and defocusing cases below (see Sec. II. B). 

B. PT-symmetric DNLS dimer 

We now turn to the PT-symmetric variant of the 
DNLS dimer, as it was recently analyzed in Ref. (28j 
(see also Ref. [22| for the experimental investigation of 
Ref. d^l). This setting is described by the system 

iu\ = —ku 2 — |«i| 2 iti — vyu\, 

iii 2 = -kui - \u 2 \ 2 u 2 + iju 2 , (5) 

which incorporates matched linear loss and gain of 
strength 7 acting on the components u\ and u 2 . The 
stationary equations have the form 

Ea = kb + \a\ 2 a + i'ya, 

Eb = ka + \b\ 2 b - i*yb. (6) 

Using the same polar decomposition as before, we find 
that the amplitudes and phase difference of the symmet- 
ric solutions satisfy 



where overdots denote derivatives with respect to the 
evolution variable t (which denotes the propagation dis- 



A 2 = B 2 = E± \/k 2 -7 2 , (7) 
sin (9 = -7/fc, (8) 
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where the signs in (JT)) correspond to the first (+) and 
second (— ) branches. These two solutions exist only up 
to 7 = k, while there is no such limit in the Hcrmitian 
case 7 = 0. The latter is the critical threshold for the 
PT symmetry breaking of the underlying linear problem, 
whose eigenvalues are E = ±\/k 2 — 7 2 . Beyond this crit- 
ical value the eigenvalues become imaginary. Addition- 
ally, as is shown in Ref. (28|, the stability eigenvalues of 
this PT-symmetric dimer are 



±2i 



2(k 2 



I 2 ) 



Ey/k 2 



for the (— ) branch and 



±2i 



2(k 2 - 7 2 ) + Eyfk 2 - 7 5 



1/2 



1/2 



system, which only allows stationary solutions with real 
E [so that | exp(iEt)\ 2 = 1]. When this symmetry is vio- 
lated, the solutions may satisfy the stationary equations 
(JBJ), but are only ghost states of the original dynamical 
system because they do not satisfy ([5]). Thus, at best one 
expects that the dynamics may stay close to the dynam- 
ics of these ghost states, especially during the evolution 
of the symmetry-breaking instability. We return to this 
topic later. 

To identify these stationary solutions, we introduce po- 
lar coordinates E — > E exp(icj)E) and get 



EA cos(4> E ) 
EB cos(<f> E ) 

Esm(cj)E) 
Ecos(4>e) 



kB cos 9 + A 3 , 
kA cos 9 + B 3 , 



k— sin 9 
A 



■7, 



A 



— k — cos 9 — 7. 

B ' 



(9) 



(10) 



To derive the asymmetric (A ^ B) solutions, we rewrite 
these equations as 



for the (+) branch. 

For the case 7 = (see Sec. II. A), these eigenvalues 
describe the critical value of E, where the two branches 
become unstable, and the pitchfork bifurcation leading 
to the asymmetric states emerges. Specifically, if k > 
(focusing nonlinearity) , the first branch corresponding to 
the symmetric solutions becomes unstable at E = 2k; if 
k < (defocusing nonlinearity) , the destabilization arises 
when E = —2k, and this happens for the second branch 
corresponding to the antisymmetric branch 9 = n. 

Similar stability conclusions occur for 7 7^ 0, where 
only the (— ) branch becomes unstable, but now also for 
7 2 > k 2 — E 2 /4. This suggests that there are two possi- Applying the identity sin 2 \ + cos2 X = 1 to pip, we ob- 



AB/k, 
2 7 AB 



cost/ 
sin 9 

Ecos(4> E ) = A z + B 
Esm(<f> B ) 



k A 2 + B 2 ' 



A 2 - B 2 



A 2 + B 2 



(11) 
(12) 



bilitics. If the propagation constant E (and the coupling 
strength k) are such that E 2 < k 2 /4, then the instability 
is induced by the increase of the PT-symmctry param- 
eter 7 at the critical point. However, if E 2 > 4k 2 , the 
instability has "already" taken place due to the presence 
of nonlinearity, and the (— ) branch is unstable even in 
the 7 = limit. In the latter case the presence of the 
gain-loss aspect only enhances the instability. 

What has become of the pitchfork bifurcation picture 
explored earlier? We can see that the same instability 
is present here (at least if E 2 > 4k 2 ). However, analogs 
of the symmetry-broken states past the critical point, 
i.e., stationary asymmetric states (emerging after the in- 
stability of the symmetric states), cannot be identified. 
Inevitably, the question of their fate arises. This type of 
question was initially raised in Ref. [l9j (where a leaky 
quantum dimer, with loss only, was considered) and, past 
a critical point, states with a complex (instead of real) 
"eigenvalue" E were identified. In the PT-symmetric 
context, a similar idea was put forth in Ref. [13] for a dou- 
ble well consisting of two delta functions. In our paper 
we unify the approaches of Refs. |lj| |34| by computing 
the "ghost states" (as we characterize them) that emerge 
from the symmetry-breaking bifurcation. 

Before presenting the computation of the ghost states, 
we comment that these states with complex (nonlinear) 
eigenvalue E are no longer true solutions of the original 
system ([5]). This is because of the U(l) invariance of the 



tain a condition for the solution amplitudes, and the same 
identity applied to (|12l) yields the parameter E. These 
solutions exist only for 7 2 > k 2 — E 2 /4. If E 2 > 4k 2 , 
they exist for all values of 7 (i.e., they have bifurcated 
"already" due to the nonlinearity). Also, these solutions 
terminate as 9 — > —tt/2, <Pe — > 7r/2. In turn, this implies 
that in this limit both B and A vanish, with the ratio be- 
tween them having the limit B/A — > (7 ± \J^ 2 — k 2 )/k. 
Thus, we have identified the disappearance point 7 2 = 
E 2 + k 2 of these symmetry-broken solutions. 

Finally, we discuss the disappearance of the two sym- 
metric states at the critical point 7 = k, which is the 
phase-transition point of the linear (and nonlinear) prob- 
lem. We have shown that at this point the symmetry- 
broken states still exist, but that now they are only ghost 
states of the steady-state problem. From the point of 
view of nonlinear theory, one may be content to find a 
saddle-center bifurcation at this point, which leads to 
the disappearance of these solutions as stationary states 
of the nonlinear problem. Yet, once again, when these 
solutions disappear (even in the normal form of such 
a bifurcation), this means that they appear somewhere 
else within the complex plane of solutions. In order to 
compare this result with the linear 'PT-symmetric case, 
where the eigenvalues collide, become complex, and con- 
tinue to exist in the complex plane, we follow Ref. [34[ 
and consider the analytic continuation of our solutions. 
In the PT-symmetric regime up to the critical point, the 
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solutions arc chosen so that u*(x) = u(—x). (This is 
a broader statement for a spatially distributed system; 
in our simpler dimer setting, we need only replace x 
by the subscript 1 and —a; by the subscript 2, or vice- 
versa.) Thus, to perform the analytic continuation, we 
use u* = us-j in j5|), which leads to ([9fl- (fl"0]) . but in the 
first pair of equations A 3 and B 3 are replaced by A 2 B 
and B 2 A. The result is 



A A 

EsinU E ) 



B _ 7± ^i 1 - k 2 
A ~ k : 

2 rji ( kl£ 



E 2 -Y T 



2 



i+e 

-it/2. 



(13) 



(14) 



Note that the pitchfork symmetry-breaking branches also 
tend to this solution, as shown above in the expression for 
B/A in the limit of termination of the branch when r ) 2 = 
E 2 + k 2 . We have made an additional subtle assumption 
here, namely, that 6 a + 9b = 0. We can obtain more 
general solutions without this assumption, but these do 
not appear to introduce new features to the problem. 

The solutions stemming from the analytic continuation 
provide a complete description of the states of the sys- 
tem. We move from symmetric/antisymmetric states to 
asymmetric ones (which may be ghost states) through a 
pitchfork bifurcation, destabilizing the symmetric (anti- 
symmetric) branch for a focusing (defocusing) nonlinear- 
ity. We terminate at the point of the linear VT phase 
transition, where the nonlinear eigenvalues of the sym- 
metric branches collide and become complex, giving rise 
to an analytic continuation of our solutions in the com- 
plex plane. All solutions terminate at j 2 = k 2 + E 2 . 



III. FROM THE NLS EQUATION WITH A 
DOUBLE- WELL POTENTIAL TO THE DIMER 

The central problem of interest here is an NLS equation 
with a double-well potential, 



iu t = Cu + iVpfU ■ 



\u\ 2 u, 



(15) 



where u t = du/dt. Here, u(x,t) is a complex field 
(which can represent the electric field envelope in op- 
tics or the macroscopic wavefunction in BECs); C = 
—(l/2)d 2 + Vi ea \(x) is a linear Schrodingcr operator con- 
taining a real, symmetric double- well potential V rea .i(x) 
(see, e.g., Refs. @, [H|); iVpj- is a purely imaginary odd 
potential with Vpj-(-x) = —Vpj-(x). 



A. Two-mode reduction 

In such settings two-mode approximations have been 
valuable tools for studying the statics, stability, and 




dynamics of the system (see rigorous justifications in 
Refs. [3(1-1381] ) . However, we do not use a Galerkin trun- 
cation to the ground and excited eigenmodes |mq, ui\ of 
the linear operator C with eigenvalues uiq, uj\ 
HI], but rather the rotated basis {ul,ur} [!, 

ul = {uo - ui) /V2, u R = (u + ui) /V2. (16) 

The subscripts L and R denote the left and right well of 
V^-oai- Note that uq and ui (the ground- and first-excited 
state) are even and odd functions of x. 

Following Ref. (39|, we approximate the solution of 
(|15l) by the Galerkin expansion u(x, t) = CL(t)v,L(x) + 
cn(t)un(x), where cl^r are unknown time-dependent co- 
efficients. Substituting this ansatz into (fTSj) and subse- 
quently projecting on ul, ur, we obtain the following 
equations for cl,r- 

icl = Ocl - wcr + i"ilcl + iil\cl\ 2 cl, 

icr = flc R -ojc L + i-fRCR + r]R}cR\ 2 CR, (17) 

where £1 = (luq + o>i)/2 and u = (wi — ujq)/2. Also, 
Il.r = J dxVvr{ x )u 2 L r, and due to the parities of Uo,i 
and Vpt it follows that 7l = — 7fi =7- To derive 
the system above we have assumed, in addition to the 
truncation itself, that the overlap integrals J dx <fr L <jr R , 
J dx<j)\4>R, and J dx<fi R (f>L are negligible in comparison 
to r/L = j dx<j)\ and tjr = J dx <fi R . This approximation 
becomes better as the distance between the wells becomes 
larger because these three integrals depend exponentially 
on the separation between the wells because of the expo- 
nential decay of the bound states Ul,r{x). (Comparisons 
of these terms with the dominant terms ?7l, t]r can be 
found in Ref. [3!|, which attests to the validity of this 
approximation.) Note that, as shown in Ref. [12j, even 
the full model with these additional terms is analytically 
tractable. However, for our present purposes, this as- 
sumption considerably simplifies the analysis and helps 
to connect with the study of the PT-symmetric DNLS 
dimer. The equality rj^ = tjr = rj mirrors the symmet- 
ric nature of V lea \(x), while 71, = — 7^ = 7, mirrors the 
antisymmetric nature of Vpf(x). 

We now consider solutions to p7|) of the form 



c L,n(t) = V 1/2 C LtR exp(-iiJt), 



(18) 



where Cl,r are unknown constant amplitudes and fi is 
the propagation constant (chemical potential in BECs). 
The prefactor r\~ x l 2 is used to rescale the term in front of 
the nonlinearity. Substituting this ansatz into (|17[) and 
using E = fi — Q, as the (nonlinear) eigenvalue, we find 
that Cl,r satisfy precisely the stationary equations for 
the DNLS "PT-symmetric dimer 



EC L = kCR + lCrfCL + ijCL, 

EC R = kC L + \CR\ 2 CR-t~fCR, 



(19) 



where we have used the notation —cu = k. 

To further justify our assumptions regarding the over- 
lap integrals, we note that for the Hamiltonian case of 
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7 = 0, the symmetry breaking is predicted to occur at 
E = — Q = 2k = — 2cj for the focusing nonlinearity case 
and at /j, — f2 = —2k = 2oj for the defocusing one. This 
implies a bifurcation at fj, = (3ujq — Wi)/2 in the focusing 
case and at /i = (3wi — wo)/2 in the defocusing case. The 
corresponding predictions for the bifurcations induced by 
nonlinearity for (|15|) arc fi ~ uq — Aq(ui — luq)/(3B — A$) 
for the focusing case, while the critical fj, for the defo- 
cusing one is fi — u} + 3B(uji — uj )/(3B — Ai) [l2j]. 
In these expressions A = J dx Uq, Ai = J dxuf, while 
B = J dxu^u\. Hence, it is clear that the approxima- 
tion for the overlap integrals made above is tantamount 
to A n = A\ = B because in that case the two expres- 
sions (the one from the dimer and the one from the ex- 
act two-mode reduction) coincide, showing that this is 
a reasonable approximation and that it trades a simple 
prediction for the critical points, which can be quantified 
just by knowing the eigenvalues of the underlying linear 
problem, for a small loss of accuracy in the result. This 
is the approximation that we use for the full problem. 



Direct simulations versus analytical 
approximat ions 



B. 



Having reduced (fT5|) to the simple dimer model of (|T9| , 
we now use the results of Sec. II to study the bifurcations 
of our PT-symmetric NLS model. To make connections 
between our work with earlier papers on double-well po- 
tentials, we use the same form of the double-well poten- 
tial that was used in Refs. [HI, H^]. It should, however, 
be evident from the exposition above that our basic phe- 
nomenology and corresponding conclusions will be quite 
general. In particular, we consider a potential of the form 
V(x) = V lea i(x) + iVpr( x )i where 

V Ical (x) = ^il 2 tI x 2 + V sech 2 (-) , 



Vpt{x) 



ex exp 



x 



(20) 



The real (even) potential V tea \{x) consists of a parabolic 
trap of strength fi tr and a localized barrier potential of 
strength Vq and width w, which has a standard double- 
well structure (that is easily achievable, e.g., in atomic 
BECs). The imaginary potential Vpj-(x) of strength e 
has a prototypical antisymmetric gain- loss profile [40| . 
Hereafter, we fix the parameters of V rca \ to be Qtr = 
0.1, Vq = 1 and w = 0.5, which yields the two lowest 
eigenvalues of the potential as luq = 0.13282 and wi = 
0.15571. This imposes a tunneling strength between the 
wells of k = -0.01145. 

There are now two possible paths to follow, both of 
which yield the symmetry-breaking pitchfork bifurcation, 
as indicated above. The first involves increasing the prop- 
agation constant /i (and hence the dimer eigenvalue E), 
thus increasing the strength of the nonlinearity. The sec- 
ond involves increasing the 'PT-symmctric potential pa- 
rameter e (and hence the corresponding dimer parame- 



ter 7). Both of these paths will yield the pitchfork bi- 
furcation when 7 2 = k 2 — E 2 /4 if the nonlinearity is 
sufficiently weak (namely, E 2 < 4k 2 ). However, only 
the latter path, which we follow below, exhibits the PT- 
symmetry- breaking phase transition when 7 = ±fc. 

In Fig. Q] we show the bifurcation diagram N = ^(7), 
where N = J^^dx \u\ 2 represents the energy in optics 
(number of atoms in BECs) . The figure compares the full 
PDE model with the PT-symmetric dimer model, where 
the real part of the chemical potential is chosen to be 
fi = 0.16. This choice of /i is justified by the requirement 
to be sufficiently close to the linear eigenvalues so that the 
two-mode picture is valid. Given the defocusing nature of 
the nonlinearity used, \i > A focusing variant of the 
problem was studied in Ref. [34j ; based on our comments 
in the previous sections, it is easy to adapt the results 
below (upon the exchange of the antisymmetric branch 
with the symmetric branch as the parent branch in the 
pitchfork bifurcation). 

There is a good agreement between the full PDE model 
and the DNLS dimer results despite the multiple approx- 
imations employed (from the original PDE to the two- 
mode system and then from the two-mode system to the 
dimer). The main features of the bifurcation diagram 
are as follows. For small values of 7 there are only two 
states with real chemical potential corresponding to the 
symmetric and antisymmetric ones for 7 = 0. These 
states are the top and the bottom one, respectively, com- 
ing from the limiting value 7 = 0. As we approach the 
PT phase transition (by increasing 7), the symmetric 
branch develops an antisymmetric imaginary part (see 
top left panel of Fig. [5]). Conversely, the antisymmetric 
branch develops a symmetric imaginary part (see bot- 
tom left panel of Fig. [2]). At the critical point the two 
waveforms arc a 7r/2 rotation of one another; thus, they 
are functionally equivalent in our U(l)-invariant setting. 
These branches collide and disappear in the saddle-center 
bifurcation, which is the nonlinear analog of the VT 
symmetry-breaking phase transition. The critical point 
for this transition is found to be at 0.01160 ± 0.00005 
from the PDE, while it is at 0.01145 in the PT-symmetric 
dimer picture. 

Before this collision occurs, the antisymmetric state 
becomes unstable (turning into the saddle of the saddle- 
center bifurcation) at the critical point 7 = 0.0083 = 
\Jk 2 — E 2 j '4. The corresponding symmetry-breaking bi- 
furcation for the PDE occurs at 7 = 0.00888. The result 
of this symmetry breaking in a purely nonlinear (non- 
TT symmetric) setting would have been the bifurcation 
of asymmetric states. Asymmetric states arise as well, 
but these are the ghost states discussed above. These 
have been computed numerically and are shown in the 
top left panel of Fig. [3] A prototypical characteristic of 
these states, shown in the bottom panel of Fig. [TJ is the 
imaginary part fij of the chemical potential [i of these 
states, which clearly demonstrates the pitchfork charac- 
ter of the bifurcation. 

The parameter fij is sclf-consistcntly computed as fol- 
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— sym 

— anti-sym-stb 

— anti-sym-unstb 
dimer-sym 

— dimer-anti-sym stb 

— dimer-anti-sym unst 
ghostPDE 

— anltic cont 

— dim-ghost 

— dim-anl-cont 
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FIG. 1: (Color online) Top panel: Bifurcation diagram 
TV = AT (7). The thin lines are the predictions of the two- 
mode analysis (the solid line denotes a linearly stable branch, 
while the dashed line depicts a linearly unstable one). The 
thick lines represent the full PDE analog (with the same 
stability designation). Bottom panel: Bifurcation diagram 
Hi — /w( 7 ). The emergence of the ghost branches is shown as 
a supercritical pitchfork bifurcation in these variables. Once 
again, a comparison to the two-mode picture is given; [ii is 
computed according to the last equation of (|12p . leading to 

m = ± 7l /(/i - + 4 7 2 - ifcVy^ - ny + 4 7 2 . 



lows. Recall that stationary solutions of (fT5|) . including 
the ghost solutions of complex [i, are governed by 



jUU = Cu + iV-pq-u + \u\ u, 
while the conjugate equation reads 

ji*u* = Cu* — iVpj-u* + |u| 2 u* 



(21) 



(22) 



Then, multiplying (f2"Tj) by u* and (|2"2")l by u, integrating, 
and subtracting the second equation from the first, we 
obtain the self-consistency condition for the imaginary 
part fij of the chemical potential: 



JZc dx M 2 



(23) 



This parameter corresponds to E sin (4>e) in the case of 
the dimer; that is, there is a direct analogy between (f2"3")l 
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FIG. 2: (Color online) Real and imaginary parts of the states 
with real (nonlinear) eigenvalues or equivalently real chemical 
potential. The left panel presents the real and imaginary parts 
of the states and their square modulus, while the right panel 
presents their corresponding eigenvalues of the linearization. 
The top panel is for the symmetric state, while the bottom is 
for the antisymmetric state. Both profiles are for the value of 
7 = 0.0100 for which the antisymmetric state is unstable. The 
two-mode prediction for the relevant instability eigenvalue is 
indicated by stars. 



with (|12j) . which becomes evident upon using the two- 
mode ansatz and the fact that 7 = J dx V-pf{x)u 2 L . In 
this variable hi one can directly recognize the pitchfork 
nature of the bifurcation at 7 = 0.00888 and the ghost 
nature of these states. In fact, only one of these is ex- 
pected to survive (the one with positive which leads 
to growth and is denoted in bold), while the one with 
negative m decays and is not be observed in direct nu- 
merical simulations. These symmetry-broken states are 
mirror images of one another, but only one has the large 
amplitude at the "right" side (for x > 0, where there 
is indeed gain) (see the top panel of Fig. . The other 
state has the large amplitude at the "wrong" side (i.e., 
for x < 0, where there is loss). 

The ghost states continue to exist past the critical 
point 7 = ±fc (as they are not subject to a bifurcation 
at that critical point). However, the other question is 
what becomes of the saddle and the center states (the 
former antisymmetric and former symmetric one) past 
the critical point of the VT phase transition. Following 
the suggestion of Ref. [34| and the fact that it is possi- 
ble to require for both of the states with real /x (existing 
prior to the VT phase transition) that u*(x) = u(—x), 
we enforce this condition for 7 > k. This provides for the 
original model an analytic continuation that is nonlocal 
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FIG. 3: (Color online) Same as Fig.O but now the two states 
have an imaginary part in their chemical potential. The top 
panel corresponds to the profile (left) and the spectral plane 
of eigenvalues (right) for the ghost state stemming from the 
bifurcation given again at 7 = 0.0100. The bottom panels, 
which closely resemble to the top ones, are for the state stem- 
ming from the analytic continuation computed for 7 = 0.0150. 



in a highly nontrivial way; the steady-state equations 
read 

fiu(x) = Cu{x) + iVpr( x ) u (x) + u 2 {x)u{—x), (24) 

and there is cross talk of the field value at x and —x 
within the nonlinear term. These states fall back on the 
original states for 7 < fc, but also provide their analytic 
continuation past this point. Interestingly, in this case, a 
new inner product is used to compute the energy (atom 
number) N. Given the substitution above, the energy 
now becomes N = dx u(x)u(—x) (see also the review 
0). Also, from Eq. wc obtain 



dx Vp T u(x)u(-3 



(25) 



The theoretical dimcr predictions for the ghost branch 
and for the analytic continuation branch are also shown 
in terms of N as a function of 7 in Fig. [T] In the 
former case of the ghost states, the normalized (based 
on the scaling transformation) norm is tantamount to 
(A 2 + B 2 )/i]l = (fi — Cl)/riL, which is a horizontal 
line shown in Fig. [T] bifurcating off of the theoretical 
point of 7 = 0.0083. On the other hand, for the lat- 
ter case of the analytic continuation states, we have 
N = f dxu(x)u(-x) = 2c L c R = 2AB/r) L = 2(fi-Q)/r] L . 
While there is a slight decay of N in the PDE for both 



states, the result is close to the theoretical prediction of 
the dimer model. A similar weak dependence of a quan- 
tity similar to A on 7 was observed in the focusing case 

Detailed examples of all branches are shown in Figs. [2] 
- [3] As described above, the symmetric branch shows 
a symmetric real part and an antisymmetric imaginary 
part, while the opposite occurs for the antisymmetric 
branch. Both branches are stable; up to the critical 
point no eigenvalues with positive real part arise. At 
that point, in the case of the antisymmetric branch, a 
pair of imaginary eigenvalues moves to the real axis (as 
shown in the bottom right panel of Fig. The vari- 
ation of the maximum real part of the eigenvalues of 
the antisymmetric branch as a function of 7 is shown 
in Fig. |U In both of the above panels the compar- 
ison with the analytical approximation given by A = 



±2iy 2(fc 2 — 7 2 ) — E\J k 2 — ■y 2 is provided, which gives 
a measure of the agreement with the two-mode discrete 
picture. 

The eigenvalues of the ghost state that emerges from 
the supercritical pitchfork bifurcation are shown in the 
top right panel of Fig. [3] Here, the solution is unstable 
too, due to a positive real eigenvalue. Importantly, one 
also sees here a shift of neutrally stable eigenvalues to the 
left-half plane, where they lead to decaying excitations. 
This spectral picture is valid for the branch of positive 
/if. It should nevertheless be recalled that this stabil- 
ity analysis is not of direct use as these states are not 
stationary, and the existence of positive fii leads to an 
amplitude that grows in time. Similar conclusions can 
be drawn for the states resulting from the analytic con- 
tinuation of the model. These states are presented in the 
bottom panel of Fig. [3] these states structurally resemble 
the ghost states, particularly in the profile of the square 
moduli. Although the model used to compute the lat- 
ter is different and involves the nonlinearity u 2 (x)u(— x), 
both states have nonzero (positive in this case) /x/, and 
both terminate at the same critical point as 7 is increased 
in the two-mode analysis. This structural proximity of 
the profiles of the top and the bottom panel of Fig. [3] can 
be interpreted as a byproduct of the predicted proximity 
for large 7 in the two-mode picture. 

Finally, we studied the dynamics of the states by means 
of direct numerical simulations. Prior to the critical 
point, we examined the dynamical evolution of the in- 
stability of the antisymmetric state, as shown in Fig. [5] 
for 7 = 0.01090. The top left panel illustrates a short 
time scale over which the mode becomes unstable. Note 
that the evolution of the antisymmetric mode between 
t = 400 and t = 500 closely emulates the growth of the 
ghost mode between t = and t = 100. Its instability 
appears to follow the growth of the corresponding ghost 
mode, which is shown in the bottom left panel of Fig. [5] 
For longer time scales, the ghost mode itself becomes un- 
stable and approaches another ghost state of the type 
that has been illustrated in Refs. [4l[ (see the bot- 
tom right panel of Fig. [SJ. The latter state involves a 



FIG. 4: (Color online) Maximum value of the real part of 
the eigenvalues as a function of 7 for the antisymmetric state. 
The instability sets in at 7 = 0.00888; for comparison, the 
corresponding eigenvalue is also evaluated from the dimer re- 
duction. 



dark soliton that has migrated to one of the wells of the 
double-well potential, accompanied by a growth of the 
amplitude and the width of the solution. Unfortunately, 
this type of state is intractable within the two-mode pic- 
ture (as the latter does not allow for the possibility of 
intra- well dynamics). In the same way, the antisymmet- 
ric state, after transiting through its corresponding (for 
the same 7) ghost state, follows the instability and fate 
of the ghost state by forming for longer times the same 
tilted dark-soliton dynamical structure (see the top right 
panel of Fig. [5]). 

Having considered the dynamical instability prior to 
the VT phase transition critical point, we now turn to 
the dynamics of various states past the latter point of 
7 = 0.01160. Above the critical point (and, in par- 
ticular, for 7 = 0.0150) we performed the evolution of 
three different initial waveforms. The first initialization 
(whose evolution is shown in the left panel of Fig. [5]) used 
the symmetric solution as obtained for a smaller value of 
7, namely for 7 = 0.0100, to examine the outcome of 
the "standard" waveforms past the "PT-phase-transition 
point which leads to the termination of their existence. 
This dynamics illustrated a phenomenology similar to 
that shown earlier in Fig. [SJ The initial growth stage 
was finally succeeded by the formation of a robust dark 
soliton on the left well of the potential. A similar evolu- 
tion shown in the right panel of the figure was observed 
for the case of the ghost state with 7 = 0.0150. We also 
performed a similar computation initialized with the re- 
sult of the analytic continuation and a similar integration 
result was obtained (not shown). We conclude that past 
the critical point of the VT phase transition, the dy- 
namics is typically attracted to dark soliton ghost states, 
which reside within the lossy well (and can thus not be 
captured by our two- mode picture). 



FIG. 5: (Color online) Contour plots showing the space-time 
evolution of the density |u(:e,£)| 2 (for 7 = 0.01090). The top 
panels show the evolution of the antisymmetric state, illus- 
trating its dynamical destabilization, and the bottom panels 
show the evolution of the ghost state. (Note the different 
time scales between the panels.) The left panels show early 
stages of the evolution (initial destabilization between t = 
and t = 400 for the antisymmetric state and then growth from 
t = 400 to t = 500). The initial stage only shows the growth 
in the case of the ghost state enabled by fii > 0. The right 
panels show the late stages of the evolution where dynamics, 
intractable within the two-mode model, develop. This evo- 
lution involves a dark-soliton state nucleated within the left 
well, similar to what was observed in Refs. [4(J. 14111. 





FIG. 6: (Color online) Left panel: Space-time evolution of 
an exact solution of the symmetric branch obtained for 7 = 
0.0100, when initialized in (|15[) for 7 = 0.0150, i.e., past the 
critical point of the transition, where the symmetric branch 
no longer exists. Right panel: Evolution of the exact ghost 
state obtained for 7 = 0.0150. 



IV. CONCLUSIONS AND OUTLOOK 

In this paper we have used analytical and numeri- 
cal methods to revisit the problem of a PT-symmetric 
double-well potential. We used a particular form for 
the symmetric real part and the antisymmetric (VT- 
symmetric) imaginary part of the potential, but we indi- 
cated that the specific choice of the potential is inconse- 
quential. In fact, given the double-well structure and a 
suitable choice of parameters (such as the chemical po- 
tential, i.e., the strength of the nonlinearity) , one can 
reduce the problem to a two-mode regime. In that set- 
ting a reduction can be made explicit (and its assump- 
tions/transformations can be suitably formulated) from 
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the original PDE problem (the nonlinear Schrodingcr 
model with a "PT-symmetric double-well potential) to 
the simpler and analytically tractable dimer setting. 

In the latter, we have shown that all characteristics 
of the system become transparent and we have identi- 
fied the symmetric and antisymmetric states and their 
nonlinear continuation. We have quantified the pitchfork 
bifurcation, formerly leading to asymmetric states due 
to nonlinearity, and presently leading to ghost states due 
to PT-symmetry. We have explicitly obtained and mon- 
itored the daughter states of this bifurcation and their 
role in the dynamics, even though they are not exact so- 
lutions. We computed the PT phase transition (or sim- 
ply the saddle-center bifurcation) of the two fundamental 
branches. Finally, the analytic continuation of the sys- 
tem past this critical point, with its unusual "nonlocal" 
(at least at the PDE level) cross talk of x and — x was 
formulated. In all of the above, the analytics performed 
in the framework of the two-mode approximation remain 
invariably a reasonable approximation to the full original 
system and the numerical computations therein. 

The above analysis gives a fundamental picture for the 
combined existence and interplay of nonlinearity, effec- 
tive discreteness (through the double-well potential) and 
PT-symmetry. However, there are numerous directions 
for interesting generalizations of this work. One of these 
would be to consider the more complicated "oligomer" 
configurations of Ref. [28j ]; namely, the trimers and 
quadrimers (for the latter, see also Ref. [42j and even the 



plaquette two-dimensonal building blocks of Ref. [43j). 
Another aspect that has emerged both here and in the 
context of Ref. [4l[ is to investigate the dynamics of 
the asymmetric states in which a dark soliton emerges 
and localizes within the lossy region. These are ghost 
states of the full PDE model; as we indicated, they can- 
not be captured by the two-mode approximation. These 
ghost states deserve a systematic investigation and clas- 
sification. Generalizing to higher-dimensional settings, 
and introducing topologically charged states such as vor- 
tices [13, EH would also be worth exploring. These topics 
are presently under study and results will be reported in 
future publications. 
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Note Added: During the finalization stage of the 
present work, the authors became aware of the recently 
posted manuscript Ref. |43| that addresses (albeit from 
a slightly different perspective) the static problem asso- 
ciated with the discussion of section II. B, namely the 
PT-symmetric dimer. 
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